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Abstract 

A shell-type model of an inviscid fluid, previously considered in 
the literature, is investigated in absence of external force. Energy dis- 
sipation of positive solutions is proved and decay of energy like 
is established. Self-similar decaying positive solutions are introduced 
and proved to exist and classified. Coalescence and blow-up are ob- 
tained as a consequence, in the class of arbitrary sign solutions. 

1 Introduction 

The following system of differential equations 

Xoit) = Vt > 

= K.iXl_,{t) - fc„X„(t)X„+i(t) Vt > 0, Vn > 1 



where fc„ = 2", has been introduced as a simplified model of 3D Euler evolu- 
tion in order to investigate a number of properties which are out of reach at 
present for more realistic models of fluid dynamics. Let us mention in partic- 
ular the works f8] [11], [12], [5], [1], [6] and references therein, devoted to this 
model and variants of it. This model differs from other dyadic or shell mod- 
els, the analysis of which is more difficult and less explicit, see among many 
others [3] , [ID] , [2] , [Z] , [1] ; system ([T]) has a basic monotonicity property that 
makes it more tractable. 

Among the many interesting properties proved in the above mentioned 
works, let us recall: i) the dissipation of energy, in spite of the fact that 



1 



formally the equation is conservative; ii) the blow-up of solutions in certain 
topologies, although the same solutions are global in larger spaces; iii) the 
relation with Kolmogorov (K41) scaling law. Our aim is to add some contri- 
bution to the understanding of these problems. The basic difference between 
the previous works dealing with energy dissipation and the present one is 
that a non-zero force was imposed until now, while we investigate the case 
without force, which contains a number of new phenomena. 

About energy dissipation, the results known until now have the following 
form. A constant positive force / > is added to the first mode 

and it is proved that there exists a unique fixed point, having Kolmogorov 
scaling, which attracts exponentially all positive solutions (solutions with 
non negative components), in the topology P. This implies that energy is 
dissipated. Positive solutions are better motivated in comparison with fluid 
dynamic equations (see [5]). 

Although the case / > and the fixed point are very interesting, it is also 
of interest to analyze the free [unforced) dynamic, namely system ([T]) without 
any forcing term. Physically, if we accept that a dyadic model like ([1]) may 
describe something of turbulence, the unforced case would correspond to free 
decaying turbulence, a widely observed phenomenon, see [9| and references 
therein. The results of energy dissipation from the previous papers do not 
extend to this case, they really require / > and make use of the non-trivial 
fixed point in the computations. For simplicity, one would conjecture the 
case / = to be similar (exponential decay to zero), maybe with a different 
proof, but it turns out this is not the case. 

For the unforced case, namely system ([T]), we first prove the following 
result: all finite energy positive solutions have energy which decays to zero 
as t — s> oo. Nevertheless, for sufficiently regular initial conditions, energy is 
conserved for small times, as was shown for instance in [8], [11] and |14j . 

Even a large class of solutions starting with infinite energy immediately 
enter namely they immediately get finite energy, and then continue their 
process of dissipation to zero. 

However, the decay of energy to zero is no more exponential. By means 
of a scaling argument we prove an upper bound on the decay of the energy 
of the order t^^. Then we prove a weak form of lower bound of the same 
order, which includes in particular the property X„ (t) dt = oo for all n 
greater than some hq. Thus exponential decay is ruled out. 
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Then we investigate positive self-similar solutions, of the form 



a. 



•n 



(2) 



t>to. 



We prove the existence of such solutions in the space The proof is 
highly non-trivial. The energy of these solutions decays exactly as and 
we conjecture that this should be the case for all positive solutions. 

The existence of such self-similar solutions relates also to the problem of 
blow-up and coalescence, but in a rather controversial way from the view- 
point of the physical interest. If we decide that only solutions with positive 
components have physical meaning, our self-similar solutions do not con- 
tribute to the problem of blow-up. If on the contrary we consider system ([T]) 
as a nonlinear model to be understood for any kind of initial condition, we 
have proved (by a simple inversion of time) that there exist solutions defined 
on the time interval (— cx),to) of the form 



with (a„) G l"^. The components X„ {t) are all negative. These solutions blow 
up in finite time, irrespective of the topology, in the sense that all components 
blow up in finite time. This result is much stronger than the blow-up results 
of the previous literature on dyadic models. 

The existence of positive self-similar solutions also implies coalescence: 
we prove that at every point of a self-similar solution there is the coalescence 
of at least another solution (which cannot be positive). 

Let us finally mention other properties of the self-similar solutions we 
construct and some open problems. 

Given tQ, we prove that there is a unique solution of the form ([2]) in 
with strictly positive components. We also prove that the components a„ 
decay as 

so, in a sense, these are 'Kolmogorov type decaying solutions' [IB]. There are 
two degrees of freedom, however, in these solutions. One is the given value 
of to. The other is that we could take 



Xn{t) 



a. 



•n 



t-to 



t<to 



ai = ... = a. 



-no — 
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and prove that there exists a unique (for given tg) self-similar solution with 
these first components equal to zero and all the others strictly positive. In 
this case, closer inspection gives us 




n/3 



where 



no — 




R e (4/5, 1) 



(numerical experiments give us i? ~ 0.885765931). 

Thus the picture is that in travel a family of self-similar solutions 
depending on the continuous parameter to and the discrete parameter riQ. 

What happens to all other solutions? Do they approach this set of special 
solutions? In this case, does a generic solution select one particular self- 
similar solution and get closer and closer to it, or does it slowly shift from one 
self-similar solution to the other? A sufficiently fast convergence to one self- 
similar solution would imply that Kolmogorov K41 scaling holds true for this 
simple dyadic model. But a slow convergence or a shift along different self- 
similar solutions could produce small deviations from Kolmogorov scaling. 
Further research is necessary to clarify these issues. 

2 Concepts of solution, their existence, posi- 
tivity and energy inequality 

To start, let us state rather general existence theorems. Their proofs are 
classical, compared to the previous literature, but at least the case of infinite 
initial energy is new, so we give the details in the Appendix. We denote by 
H the space namely the space of all square integrable sequences of real 
numbers. It is a Hilbert space with the obvious inner product {u,v)fj = 
Y^n=i'^nVn, whcrc u = (Mn)„6N' ^ = (^«)nGN' u,v G H. The Corresponding 
norm in H will be denoted by We denote the space of all sequences of 
real numbers by M^, and its subset of all non- negative sequences by M^. 

In the sequel we shall often use the term "energy" for the quantity 
for an element X & H oi also X G (it may have infinite energy). 

Definition 1 Given X^ G M^, we call componentwise solution of system ([1]) 
with initial condition X° any sequence X = (X„ (■))neN ^/ continuously dif- 
ferentiable functions X„ (■) : [0, oo) — > M such that X„ (0) = X° for alln & N 
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and all equations in system ([T]) are satisfied. If X (t) G for all t > 0, we 
call it a positive componentwise solution. 

If X is a componentwise solution, from system ([1]) we have 

ft 

Xnit) = xOe-'^-^o ^ / e-^"^>"+i('^)'^"/t„_iX^_^ (s) ds. (3) 

Jo 

This identity will be used quite often. 

Definition 2 We say that a componentwise solution X has finite energy for 
positive times if X (t) G H for all t > 0. If also X° G H, we call X a finite 
energy solution. 

Theorem 3 Given X^ G M^, any componentwise solution of system ([T]) with 
initial condition X^ is positive. At least one such solution exists. Moreover, 
any such solution has the following properties: 

i) for every n > 1 and t > we have 
and hence 

n n 

j:x^(t)<Y.{x^r (4) 

i=i j=i 

a) if X^ > for some n > 1, then X^it) > for all m > n and allt>0. 

Theorem 4 For every X^ G H, there exists at least one finite energy solu- 
tion of system ([T]) with initial condition X^ , with the property 

\X (t) 1^ < \X{s) \h for allO<s<t. (5) 

Moreover, if X^ G iJnlR^, then all componentwise solutions are finite energy 
and satisfy ([5]). 

See the proofs in the Appendix. 
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3 Energy dissipation 



System ([T]) is formally conservative: 




n=l 



by a simple rearrangement of the series and the condition Xq = 0. This 
rearrangement is rigorous if the solutions live in a sufficiently regular space 
(see for example [5]). Such kind of regularity may hold for small times if the 
initial condition is very regular itself (see [S]), but for sufficiently large times 
we prove that solutions dissipate energy. The intuitive mechanism is a very 
fast shift of energy from small to large n components. 

We give two results of energy dissipation for positive solutions: infinite 
initial energy becomes finite immediately; the energy of a finite energy solu- 
tion tends to zero as t ^ oo. Although the degree of infinity of the energy 
of initial conditions can be generalized, for the simplicity of statements we 
restrict ourselves to X° of class the norm ||X°||oo := sup„ \X^\ is finite. 

Theorem 5 Assume X^ G fl and let X be a positive componentwise 
solution of system with initial condition X^ . Then X has finite energy 
for positive times. 

Theorem 6 If X is a positive finite energy solution, then 



Moreover, given L > and a > 0, there exists t > depending only on L 
and a such that for all positive finite energy solutions X with \X{0)\h < L 
we have \X (t)|^ < a. 

The proof of both statements is based on the following lemma 

Lemma 7 Let X be a positive componentwise solution, with X (0) G /°°nM^ 



and let ||X(0)||oo < L. Let (j)n{t) := ELi^lW f^m > 1, let (f)o{t) = and 



Then there exist two summable sequences of positive numbers {a„}n>i and 
{sn}n>i depending only on L, such that: 



lim = 0. 



t— >oo 



Mt) ■■=EZiXm = \x{t%. 
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i) for all n > 1, for all t > and for all e G (0, 1] one has 

(t) < ea„; (6) 

ii) for all integers M > 1 one has 

Proof. For all n, let := 0„(O)VL, so that by equation (jlj), X„(t) < rrin 
for all t > 0. 

We observe that it is possible to find two summable sequences {a„}n and 
{Sn}n such that for all ri > 1: 

2"s„a^ > 2m^m„+2 ( 1 ^ 



It is enough to set s„ = 2~"/^ and a„ = (72""/^ with C a suitable constant, 
and recall that by hypothesis rUn < Ly/n. We observe that without changing 
the sequences, a fortiori for all 5 < 1, 



2"^n< > 2m>„+2 1 + ^_22., ^ • (8) 

Part 1. Given the two sequences, we now show that the upper bounds (E]) 
hold. 

Let h = e^'^Sn- Since is nonincreasing, for all s G [0, h], 

(f)n{t + h)- (j)n-l{t) < Mt + S)- <Pn-l{t + s) = Xl{t + s), 

hence, if for some s inside the interval, X'^{t + s) < £:a„ we are done. 

On the other hand, let us suppose that X'^{t + s) > can for all < s < /i. 
One has 

(j)n{t + h)- = Xl{t) + / —Mt + s)ds. 



ds 



We need to prove that 

d 



, -^n{t + s)ds < ean - X^{t). 
ds 
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For sake of notation simplicity, let a = eankn and let A = kn+iTnn+2- Then, 
by Theorem [3]- 

'n{t + s)ds= I -knXl(t + s)Xn+iit + s)ds < -a Xn+i(t + s)ds. 



ds 

We need a lower bound for Xn+i- In the interval [t; t + h] we know that 



Xn+l — knX'^ — kn+lXn+lXn+2 > « " AX„+i, 



whence we get 



Xn+i{t + s)> e-^^X„+i(t) + ^ (1 - e-'^) > ^ (1 - e-^') 



By substituting into the integral above one gets 

I ^Mt+s)ds < (1 - e-^^) ds = -y,{e-''-l+Xh) < "^^YT^ 

Substituting next h, a and A and recalling the condition ([8]), we finally get 



-^(1 + 2m„+2 ( 1 + 



This concludes the first part. 

Part 2. Let M > 1 and define the sequence {t„}„>A/-i by tA/-i = 
and tn = e~'^ Ylt=M ^k, for n > M. We substitute t = t„_i inside the 
inequality IQ: 

4>n{tn) - 4>n~l{tn-l) = 4>n{tn-l + - <^„_l(t„_l) < £a„. 

Adding the above inequalities for n from M to any number N > M, one has 

N N 

0Af(tAf) - 0M-l(O) = ^ [0n(in) " 0n-l(in-l)] < ^ ^ On- 

n=M n=M 

Monotonicity of 0Ar yields that 

AT 



(£ ^^Sfc) < (pNitN) < 0M-l(O) + £ On, 



k=M n=M 



hence letting go to infinity we get ([7]). The proof of the lemma is complete. 
■ 

We are now ready to prove the above theorems. Consider the assumptions 
of Theorem [51 By letting e = 1 in the second part of Lemma U\ one has that 
X{t) G H for t > Yl'k=M ^k- Letting M go to infinity we get the thesis of 
Theorem O 

As to Theorem [HI let us prove that for all a > there exists some t > 
such that 0oo(t) < a. Since ||X(0)||oo < \X{Q)\h < L, we let M = 1 in 
the second part of Lemma [3 we only need to choose e in such a way that 
e Yl'^=i < «• One gets i = e"^ YlkLi ^k- The proof is complete. 

4 Bounds on the decay of energy as t ^ 00 

In this section we prove a bound from above and another from below, for 
the decay of energy as t — > cxd, which essentially say that solutions decay as 
t~^. The results are restricted to positive componentwise solutions. The first 
result is due to a scaling argument based on the fact that the nonlinearity is 
homogeneous of degree two. 

Theorem 8 Let X be a positive componentwise solution, with X (0) G fl 
M^. Then there exists C > such that 

fort > 1. 

Proof. First, by Theorem \5[ the solution has finite energy for positive 
times. Hence, by Theorem [Q there is a time to such that \X (to)|/f ^ 1- It is 
thus not restrictive to prove the theorem in the case jX''!^ = 1. 

First a general fact: from Theorem [6] we know that there exists i > 
such that for all initial conditions X° with = 1 we have 

\Xii)\l<\- 

Let us start the proof, for a solution X such that jX*^!^ = 1. Let = 
\X Consider the rescaled if-valued function Y = (F„)^gj^ defined as 

Y{t) := a^^X {a^H + t) . 
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We have 

Yn{t) = afK-iXl_^ {a^H + afknXn {a^H + i) X^+i {a^H + t) 

= (t) - KYn (t) Yn+l (t) . 

This means that F is a finite energy solution of system ([T]). Moreover, 
|F(0)|5^ = 1. Hence 



ttn := 



1 



\yit)\H<-.- 



In terms of X we have 



X[-+i 



2 2 
^1 ^2 * 



H 



By induction we can prove that 

i 



X 



t 



— + ■■■ + — +t 



for every k>l. Since < ^, we have 



2 2 



t t - 2t 

+ ••• + — + t< . 

ai ■ a2 ■ ■ ■ ■ ■ dk Oil ai ■ a2 ■ ■ ■ ■ ■ Oik 

Recall that energy inequality holds for all positive finite energy solutions, see 
Theorem O Hence for all t > 



2t 



we have 



If we restrict to — 

ai-a2-----afc 

Hence 



< t < 



2t 



ai-a2-...-afc-Q!fe 

|A-WP„<f 



— we also have af ■ ■ . . . ■ o:]:,^ < 



2t 



for all ^ < t < - 

The proof is complete. ■ 



This implies the claim of the theorem. 



Theorem 9 Let X be a positive componentwise solution, with X (0) G fl 
M^. Let no + 1 be the minimum integer with the property ^^^+1 > 0. We 
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know from Theorem\^ that > for all n > Uq and t > 0. Then, for 

some constant C > 0, for every n > Uq and t > 1 we have 

II X„+i(.)rf. > logt + k-' log ■ 

Thus, in particular, for every n > hq, 

/oo 
Xn+i{s)ds = OO 

and 

limsupt ■ Xn+i(t) > k~^. 

t^oo 

Proof. From identity we have 

By the upper bound of the previous theorem, there exists a constant C > 
such that 

hence 



namely 



-^n I Xn+l{s)ds < log (^Y^i)T 
rt log (^)+ logt 

/ X^^,{s)ds> ^ . 



This implies the claims of the theorem. The proof is complete. 



5 Self-similar solutions and related facts 

We call self-similar any solution X of the form X„(t) = a„ ■ <f{t). It is easy 
to verify that self-similar solutions satisfying the equations have the form 

Xn{t) = 7^, t>to, (9) 
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with to < 0. We are interested in finite energy self-similar solutions, hence 
we require also 



E 



< oo. 



n=l 

In the next section we prove the following result. 

Theorem 10 Given tg < 0, there exists a unique finite energy self-similar 
solution with Oi ^ 0. In general, given to < and Uq > 0, there exists a 
unique finite energy self-similar solution with 

ai = ... = = 0, ttno+l 7^ 

(where the first conditions are meaningful only for Uq > 0). In addition, the 
coefficients an have the property 

_i /c. — CnQ. 

n^oo U -^/•^ 

Thus we see that Kolmogorov scaling law [T3j (so called K41) appears 
in these special solutions, phenomenologically associated to decaying turbu- 
lence. But it is for us a very difficult open problem to understand whether 
all other solutions approach the self-similar ones and in which sense, so we 
cannot say how general this scaling property should be considered. 

The existence of finite energy self-similar solutions is of conceptual in- 
terest in itself, in comparison with analogous investigations for Euler and 
Navier-Stokes equations. Also here, in this very simple context, the proof 
is highly non trivial. Apart from its intrinsic interest, the existence of such 
solutions has a number of implications. 

First, they realize perfectly the decay t^^, coherently with the previous 
section (where the lower bound was more vague). We conjecture that the set 
of all finite energy self-similar solutions, set depending on to G M and no > 0, 
attracts all other solutions in a suitable sense. If this is the case, the decay 
t~^ would be more strictly the true one for all solutions. 

A second implication is the existence of solutions that blow-up backward 
in time. This is of interest for two reasons. To explain them let us first clarify 
what happens to solutions when we reverse time. 

We may consider system ([T]) for negative times, t < 0, and give a def- 
inition of componentwise solution exactly as for t > 0. All definitions and 
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theorems can be rewritten for negative times. One also has the following cor- 
respondence between the forward and backward problem: let X = {X (t))^^Q 
be a componentwise solution of system ([1]), for t > 0, as usual. Then 
Y = {Y (t)),<o defined as 

Y (t) = -X i-t) (10) 

is a componentwise solution for t < 0. Indeed, by the purely quadratic nature 
of the equation, 

Y^it) = = K-lX'n-li-t) - KXn{-t)X^n+l{-t) 

= kn^iY^_^{t) - knYn{t)Yn+l{t). 

Thus, any positive solution over [0, oo) gives rise to a negative solution over 
(— cx3,0], and vice versa. 

Theorem [in] ensures that there exists a self-similar solution := 
with to < and a„ > for all n > tlq. It is easy to check that Xn{t) is a 
componentwise solution on the open interval (to, +C)o). The energy is finite 
for all t > to and limj^oo |-^n(t)| = whereas \im^^^+ Xn(t) = +oo for all 
n> riQ. With time inversion ( ITOl) it is possible to define Y{t) := —X{—t) for 
all n and t G {—oo, —to)- Y{t) is a componentwise solution of ([1]) on the open 
interval (— oo, —to) and limt^_oo |^n(t)| = whereas \im^^_^- \Yn(t)\ = +oo. 
This means that Yn blows-up in finite time. Notice that every component 
blow-up, not only some norm of the solution. 

A consequence of this is a coalescence property of self-similar solutions. 
Theorem H] applied with initial condition 1^(0), states that there exists a 
finite energy solution Y{t), defined for all t > 0, with bounded energy and 
initial condition 1^(0) = Y{0). Moreover, the condition ([5]) of the theorem 
ensure that the energy of Y is nonincreasing in t, on the contrary \Y\h is 
increasing in time. This means that on the system ([T]) there is not uniqueness 
of solutions in M^. This result is of interest in itself, but we stress that it 
holds in the enlarged class of non-necessarily positive solutions. 

Let us invert the time again with X]sf{t) := —YN{—t). X and X are two 
componentwise solutions on the interval (to,0], with decreasing in time and 
nondecreasing in time energy respectively. They are different on (to, 0) and 
coincide at time t = 0. The role of time t = in the previous argument can 
be replaced by any time ti > to. Thus we have 

Corollary 11 If X is a self-similar solution of the form (Q, all its values 
are coalescence points, in the sense that for all ti > to there exists a finite 
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energy solution X^^ , defined fort G (— oo,ti], such that X*^^ (ti) = X (ti) and 
X*i(t)^X(t) on (to,ti). 

Finally, we state the blow-up for negative self-similar finite energy solu- 
tions. 

Corollary 12 Given to > and uq > 0, there exists a unique negative finite 
energy self-similar solution, defined on [0, to) with 

ai = ... = a„Q = 0, a„g_|_i ^ 0. 

This solution blows-up at time to. 

Proof. Apply the inversion flTUl) to the solution given by Theorem [TUl ■ 



6 Existence and uniqueness of self-similar so- 
lutions 

In order to prove Theorem [10] it is best to restate the problem in terms of 
properties for sequences of positive numbers. 

If a positive componentwise solution is of the form (Q, then 

(3 ci^ a, Oj 

— Xnit) = kn-lX^_^(t) — knXnit)Xn+l(t) = kn-l' " 



(t-to)^ — -^n-iv^. -\t-to) 

SO the sequence {a„} must satisfy 

for all n. It is indeed possible for the first terms ai, 02, . . . , ctno to be zero, 
but by induction, if a„(,+i > then the subsequent coefficients must satisfy 

a„+i = 2-" + ^ > 0, n>no + l (11) 
To simplify matters and take into account no, let 
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so that the condition becomes = 0, Oi > and 

ttn+i = 2 + 4^, n > 1 (12) 

a„ 

Since, given ttno+i or ai, this recurrence uniquely defines the sequence, in 
order to prove Theorem [10] we have to show that there exists a unique positive 
number ai such that the sequence {an}n G H. We will prove the following. 

Theorem 13 There exists a unique real number 7 such that the sequence 
{a„}„ is in H iff ai =7 ( equivalently , iff = ai = ■ ■ ■ = = and 
ano+i = 2-"«-i7;. 

Moreover, let (3 = 2^^^^ . One can find R > and a strictly decreasing 
bijective function h : (0; R] — > [0; 00) such that 

z; 7 = hi(3'R), 

a) if {a„}„ G H , then an = 2~'^h(P'^^^~'^^^ R) for all n > uq. In particular 
an ~ , for n 00, with Cn,, = P'^^'/R. 

Numerical computations give 7 ^ 0.917576296. 

It is clear that Theorem [10] follows immediately from Theorem [131 The 
proof of the latter requires some work and will follow from Theorem [T7] below. 
Three technical lemmas will be needed. 

We start by looking for another sequence {dk}k>-i (not depending on 
a„), satisfying the peculiar relation below: 

rd_i = -i 

I Y,d^dk-^i^-/3'''^'')=2dk0■ 

Lemma 14 Equations f[T3l) uniquely define a sequence of real numbers {dk}k>-i 
such that dk > for k > 0. The power series YlT=o^kx'^ ^ positive con- 
vergence radius R > 0. If we let h{x) := —YlT=-i^kx'' , the function h is 
defined on the interval (0; R) where it is analytic, nonnegative and strictly 
decreasing, with h{0~^) = +00 and h{R~) = (so that it can be continuously 
extended on {0; R]). 



k> -1 



(13) 
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Proof. The second one of 0131) uniquely defines d^^i as a function of the 
previous terms. Truly, the coefficient of d^+i is (i_i(8 — /^^'^"^ — which 
is always nonzero. 

With some algebraic manipulations, the above system of equations can 
be rewritten as: 

= -1 

rfo = (2/?-^ - /3 - ~ 0.8155665 (14) 
^dk+i = \J2i=o^k,ididk-i, k>0, 



where 



f 1 + /^^ - 

1 - /?7 _ fjU 

(l-/3^^+^)(l + /j^^+^) 

1 - /32fc+7 _ /54fe+ll 
X _ p'ik-2i+9 _ p2k+2i+9 



k = 0, i = 
/c > 0, i = 0,k 



<i < k 



1 _ f]2k+7 _ pik+U 

All ak,i are positive and so (ifc > for k > 0. We notice moreover that 
lim/j^oo Cfcfc.i = 1 uniformly in i (exponentially in k, see the Appendix). 

The recursion (IT^ is similar to the classical Catalan sequence, but since 
we have only some asymptotic control on the coefficients a^^i, and since the 
behaviour of {dk}k strongly depends on the ffist values (even its convergence 
radius does), we need some explicit investigation of its properties. 

For k > 0, let d'l^^^ := | Yl^=oi'^k,i — ^)didk_i and let d'^ = do. 

Let g{x) := Yl'^=o ^kx'' so that h{x) = —g{x). Let g{x) := Yl'k'=o ^'k^^ ■ 
By the third one of (|T4|) it follows that 

oo ^ oo A: 

dk+ix^^^ ^ ^2 ^ ^ '^k,idix'dk-iX^~^ 

k=n k=n i=0 

If n is such that \ak,i — 1| < e, we have 

n ^ oo fc ^ ^ 

9{x) - ^ d'f^x^ < dix'dk-ix''-' = -^xg^{x) 

k=0 k=0 i=0 

Writing also the corresponding lower bound and letting e — > and n — >■ oo 
accordingly, yields 

g{x) -g{x) = -xg^{x). 
16 



The above formula is true for all (complex) x inside the radius of convergence 
of g (in the Appendix we show that g has a convergence radius times 
larger than g). 

For all X inside the convergence radius of g^ g{x) must be one of the roots 
of the above degree 2 polynomial, i.e. 

1 - a/1 - 2xg{x) 
g{x) = ^ 

X 

(The other root does not satisfy 5^(0) = do.) 

In the Appendix we show that the radius of convergence of g is greater 
than 1 and that inside the unit circle of the complex plane l — 2zg{z) has only 
one zero, which is in fact some real point R G (4/5, 1). It follows that R is the 
radius of convergence of g, that h is defined on {0]R) and that h{R~) = 0. 
The other properties follow from inspection of the first coefficients dk- ■ 

From now on h is intended to be extended up to R. Since h is bijective 
and has image [0, cxo) and recalling that a„ > for all n, it makes sense to 
compute h~^(dn)- 

Lemma 15 Let {a„}„>o be any positive sequence satisfying (fT2|) . not neces- 
sarily with flo = 0. For n > let A„ := h~^{an) ^ (0,-R]. The numbers A„ 
satisfy: 

Proof. Recall that h{x) = —J2T=-i^kx'' ^ 

oo oo fe+1 

fln-i = h'^i^n-i) = ^ didjXl^_i\i^_i = ^ A,^_i ^ didk-i 

i,j=~l k=-2 i=-l 

Using relation ([T5]l for /c > —1 and since /52A:+2i ^ 4 jf ^ _ _2 -^^^ get: 

00 00 k+l 

k=-l k=-2 i=-l 

00 00 

= 2 J2 dkf3''X'n-i+ Yl I3^'^^'d,dj\';^\ = -2/i(/?2A„_i)+/i(/9'A„_i)/i(/52A„„i) 

k=—l i,j=—l 



(15) 
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Substituting in equation (fT^ we find 

which is equivalent to ffTB]) . ■ 

The ratio between A„ and /?^A„_i plays a crucial role in what follows. 
Next lemma characterizes its logarithm. 

Lemma 16 Let {a„}„,>o be any positive sequence satisfying (fT2!) . not nec- 
essarily with Qq = 0. For n > 1 let := logA„ — logA„_i — log/9^. The 
sequence A^ can present two behaviours only: 

1. If for any n>l, A'„ = 0, then A^ = for all m > 1. 

2. //A^ 7^ for some n>l, then there exist two constants c,d> such 
that for all k > 0: 

\XU\>c{l + d)'\X'J, {-1)'X'^_,X>0 
Proof. If A^ = 0, then of course A„ = /3^A„_i, hence ([T6|) reduces to: 

Since h is injective, A^+i = /3^A„, that is X'^^^ = 0. This proves by induction 
the first part of the lemma, for k > n. If the same were not true for some 
k < n, then by the second part it would have been A^ 7^ 0. 

In particular note that if = (yielding Aq = R) and Ai = P'^R (yielding 
X'l = 0), then by the first part and recurrence ( |T2l) . h{P^R) = 02 = 2. This 
will be used in the second part. 

For proving the second part, let ip{x) := log/i(e^') and let Ax{y) : = 
h{xey)/h{x). 

The function is defined on (— cxd, logi?), where it is analytic, strictly 
decreasing, concave and invertible; ip has one oblique and one vertical asymp- 
tote: ip{x) + a; t for X I —00 and ip{x) ~ | log(log-R — x) for x t logi?, so 
that the image of ip is the whole M. Moreover ip' < —1 on all the domain. 

Ax{y) on the other hand is defined for x G (0; R) and y G (— cc; logi? — 
logx) and this last interval always contains a neighbourhood of 0. A^(-) is 
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analytic and decreasing on all its domain, moreover we have: 

A,(0) = 1 

>e-y>i, y < 

AM < < 1, ?/ > 

in fact, 

^x{y) = h{xey)/h{x) = exp{tjj{\ogx + y) - ^(logx)) = exp{i;'{^)y) 
and ip'i^^) < —1. By virtue of (fT6l) we can write: 



n+l) 



n-l 



1 - 



KXn) 



+ 



h{(3^X^_,)h{l3^Xn_,) 

hiP^XnMXn) 



hi(3^Xn-l) hiP^Xr. 

j^^^^l^ ^) A/^-a.(-a;) (1 - a,„(-a;)) + a,.,„(-a'ja,„(-a;) 

'1 1 
2 



Hf3''X^-i 



a^.,„(-a;)a,„(-a;) + 



Let 



1 



1 

< - 

- 2 



KP^Xn-i] 

We need to prove that Kn is always nonnegative, and bounded below by a 
positive constant for n large enough. 

Since A„_i < R, h{p^Xn-i) > h{p^R) = 2. So > for all n > 1. For 
n > 3, a„_i > 2 by (fT2l) . hence A„_i < /?^-R, so: 



A-„ > i 



1 



1 
2 



1 



--■.K>0. 



Going back to the previous expression, if A^ > 0, since both coefficients are 
positive, 

A/3.A„(AUi) = 2ir„A^.,„(-A:jA,„(-A:,) + (1 - 2if„)A^.,„(-A:,) 
> 2Kne'^^'- + (1 - 2Kn)e^" 
= 2Kne^^'" + (1 - 2/s:„)e^'" - e^" + e^'" - 1 + 1 
= 1 + (2is:„e^" + l)(e^" - 1) > 1 
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A posteriori we get A^^;^ < and so, after letting x := e"**", 
e-^"+i > A^2;,,^(A;+i) >2ir„x2 + (l-2ir„)x=:ir„x(x-l) + e„(x) >^^„(x), 
where 

Onix) := KnX'' + (1 - Kr.)x > x'^""- =: Onix) 

The above inequahty holds for all x > 1 and follows by comparing the two 
functions, their derivatives and by the fact that K„ < 1, in fact: 

6,^1) = 1 = ^„(1) e'^il) = ^+Kn = O'nil) 

e'^ix) = 2Kn >Kn + Kl> {Kn + K^)^^^"-^) = Ol^x) , ioi X > I 

Finally e"'^"+i > 9n{x) = e'^^+'^")^" and hence 

A;>0 ^ A;+i<-(l + ir„)A; 
On the other hand, if A'„ < 0, with analogous reasoning: 
e^K+i < Af32x„{X'n^i) < 2KnX^ + (1 - 2ir„)x = 2K„x(x - 1) + x < x = e^'" 
and hence 

A^ < ^ A^^i > — A^ 

Summing up things, for any A^ ^ 0, the sequence {A^_,_;.}fc has alternating 
signs, moreover: 

\X'n+k+2\ > (1 + -^n+fc)|A^+^| 

The second part of the lemma is now proved by induction on recalling 
that Kn>K>{]ioin>?,. ■ 

We are finally able to make the main statement. 

Theorem 17 ^4 sequence an satisfying f|TT]) with an^ = and a„p+i > can 
present only two behaviours: 

1. If for any m > 1, A^ = 0, then there exists c > such that for all 
n > uq, an = 2~"/i(/5^("~'^°)i?); in this case a„ ^ for n f oo and in 
the Urmt an = CnoP"" + 0(/33"), with = /J^^o//?. 

2. If for any m > 1, A^ 7^ 0, then either the odd or the even terms of an 
diverge more than exponentially and the power series anz"' has zero 
radius of convergence. 
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Proof. We apply directly Lemma [TBI 

In the first case, for all n > 1, = 0, hence A„ = /?^"Ao = /3^"-R, by 
which ttn+no = 2-"o-"a„ = 2-"o-"/i(/32"i?). The limit behaviour of a„ follows 
from h{x) = + 0(1) for x i 0. 

In the second case, notice that for all k >2, 

logak = log/i(Afc) = ^(logAfc) = V^(logAfc_i + log (3^ + A'J > ^(A'J 

Let n = m or n = m + l so that A^ < and hence A^_,_2fc < for all k and 
Ki+2k ~^ k ^ CO. By the properties of ip, there exists Xq such that 

X < Xq ^ 4'{x) > —X — 1. We deduce that, for k large enough: 

loga„+2. > m.+2k) > -K+2k - 1 > c(l + df^\\'^\ - 1 
so the power series 'Ylik^^z^ radius of convergence and the same for 

Proof of Theorem 1131 By Theorem [T71 we are in the first case iff A'^^ = 
that is if ai = h{(3'^R). In the first case {a„} G H, since a„ = C^g/S" + 0(/5^") 
and (3 <1. In the second case of course the sequence is not in H. 



7 Appendix 

7.1 Proof of Theorem [3] 

The first claim is an obvious consequence of identity ([3]). 

Let us prove the existence statement. Given X° G M^, consider the finite 
dimensional system 

Xoit) = Vt > 

Xnit) = kn-iXl^^{t) - A;„X„(t)X„+i(t) Vt > 0, G {1, 2, . . . , N} 

X^+i(t) = Vt>0 

X„(0) =X0 VnG {l,2,...,iV} 

It has a unique solution: local existence and uniqueness comes from Cauchy 
theorem (the vector field on the right-hand-side is locally Lipschitz contin- 
uous), global existence is a consequence of the bound on maximal solutions 
derived from 

N N 



Y.^i{t) = Y.^m (17) 



n=l n=l 



21 



which is easily proved by computing ^^n=i-^ni't)- Denote by the 
unique solution. 

Since X° G M^, the solution is positive, in the sense that X^{t) > for 
every n = 1, ...,N and t > 0. This simply happens because 



Moreover, for the solution X , for every n G {1, we have 



dt 



Hence, being positive, we get ^ YTj=i {-^fY — namely 



j=i i=i 



In particular. 



< (t) < 



(18) 



J2 {XjY for alH > 0. 
i=i 

On a bounded interval [0,T], consider now the family (X^)^^^, for a 
given n > 1. The assumptions of Ascoli-Arzela theorem are satisfied for 
this family. Indeed, equi-boundedness has been proved above in f[T5]) : equi- 
uniform-continuity (equi-Lipschitz continuity, in fact) comes from the iden- 
tity 

ft 



it) -x:: is 



kn-i{X^^,{r)Y - A:„X„^(r)X„^^,(r) 



dr 



and the already proved equi-boundedness (recall always that n is fixed). 



Thus, from Ascoli-Arzela theorem, for every n there exists a sequence {N^"^} 



ken 



such that {Xn * jfceN converges uniformly to a continuous function X„. By 
a diagonal procedure, one can modify the previous extraction procedure and 
get a single sequence (Xfc)^^^ ^^ch that (X^*)^^^^ converges uniformly to 
Xn- We can thus pass to the limit, as ^ oo, in the equation 



X^^it) =X:+ fc„_i(X„^^,(r))^ - A;„X„^Hr)x5, 



dr 
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and prove that 

X„(t) =Xl+ f {K_,Xl_^,{r) - A;„X„(r)X„+i(r)) dr. 
Jo 

Thus the functions are continuously differentiable and satisfy sys- 

tem ([T]). Of course they are non negative, being the uniform hmit of non 
negative functions. Continuation from an arbitrary bounded time interval to 
alH > is classical. 

Finally, we have to prove properties (i) and (ii) of the theorem, for any 
positive componentwise solution X. As to (i), we repeat the argument used 
above for X^: for every n > 1, from system ([T]), we have 

d ^ n 

— ^ X| = 2 ^ kj^iX'^^_^Xj - kjXjXj+i = -knX^Xn+i 
j=i j=i 

hence, X being positive, we get ^ Sj=i {-^j^Y — ^^us (i) is proved. 
As to (ii), from ([3]) we have 

Xn{t) > xOe-'^-^o^^+iM*' > 0. 

For Xn+i(t) we use the inequality (again a consequence of ([3])) 

Xn+l{t)> / ^'^^"+^^'''^'^''knXl{s)ds 

Jo 

and the fact that X„(t) > for all t > 0. By induction we get the result. 
The proof of the theorem is complete. 

7.2 Proof of Theorem [4] 

We introduce the same finite dimensional system as above. We do not have 
anymore the positivity property, but we still have (fT7|) . Since the (global) 
initial energy is finite, we have 

N 

j:{x-it)r<\x% 

n=l 

for alH > and > 1. This implies again a bound on single components: 

\Xnit)\<\X'\^ 
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for every t > 0, > 1 and n = 1, N. Having this bound, we proceed as 
above and prove, on a given [0, T], the existence of a componentwise solution 
X, with X^'' Xn uniformly on [0, T] as — > oo, along a diverging sequence 
(iVfe)fceN- Fi^om 



n=l 



one easily get 



2 
H 



n=l 



first by truncating the sum up to a given value R and taking the limit in k, 
then sending R to infinity. Hence in particular X {t) E H for alH > 0. 

Finally, assume that X° E H n and let X be any componentwise 
solution. From Theorem [3] (i), it satisfies 



j:(x,(t)f<Y,{xn-<\x''\i 

for every n > 1 and t > 0. This implies that X is finite energy and satis- 
fies ([5]). The proof is complete. 



7.3 Estimates on ak^i 

Recall the definition of akX- 



OikA 



1 - - 

(l-/32^'+2)(l + /?2fc+7) 

1 - /?3fc+6 cosh((A; - 2i) log/3) 



2k+7 



jj4k+n 



k = 0, i 
k > 0, i 
<i < k 





0,A; 



We study their behaviour for k large and look for bounds. Notice that the 
denominators are all positive for A; > 0. 

Lemma 18 There exists C > such that for all k >0 and all < i < k , 

<C/32'. (19) 
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Proof. There are three cases. 
1) When A; > 1 and i — or i — k, ak,i < 1, in fact: 



1 - OikA = 



1 - /32fc+7 _ /54fc+ll 
1 _L /^2fe+7 



1 - /32fe+7 _ /34fe+ll 



which is positive. The last factor is decreasing in k and is equal to 4/3 
when A; = 1, so for A; > 1 

Q<(3^-(3'< (3-'\l - < 4/3(/52 - (5') =: C 
2) When A; > 2, < i < /c, we have a^^i > 1, since, recaUing that cosh(a;) < 

-/33'=+6cosh((A; - 2i) log/5) + + /J^'^+i^ 



> 
> 



1 - /52fc+7 _ ^4fe+ll 
_ j^Zk+& j^-\k-'2i\ _|_ p2k+7 _j_ ^4fc+ll 

1 - /?2fc+7 _ /J4fc+ll 
_p2k+8 _j_ p2k+7 _j_ ^4fe+ll 
1 - _ /34fc+ll 

2k+7 l-P + P'"^' 



^ I — p2k+7 _ p4k+ll 

The right factor is decreasing in k, hence for A; > 2 

/3-2'^K,-l)>/3^-/?«>0 

Prom cosh(x) > 1 it follows that 

_p3k+6 _j_ p2k+7 _j_ p4:k+ll 
Oik,i — 1 < ]^ _ p2k+7 _ ^4fe+ll 

^ ,2k-,7 1 - P'-' + P''^' 
^ 1 _ /32A;+7 _ pAk+U 

The right factor is increasing in A; for A; > 2 and it tends to 1, so 
0<P'-P^< p-'\ak, -1)<P'<C 

The third case, namely A; = is included by adjusting C. ■ 
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7.4 Estimates on the convergence radius of g and g 

By equation f|T5l) . in order to find the radius of convergence of g, we need to 
know tlie zeros of g{z) := 1 — 2zg{z). 

Theorem 19 The complex function g has only one root inside the unit circle 
centered in the origin. This root is a real number R G (4/5, 1). 

The proof is based on Rouche's theorem for holomorphic functions. 

We spht g in the sum of two terms: g = gA + ds, where ^a(^) = 1 — 
2dQZ — 2d[z^ and gsiz) = —2z Yl'^=2 ^'k^''- Rouche's theorem tells us that if 
IoaI > Igsl on the contour of some compact set of C, then cja and g have the 
same number of zeros inside the compact. 

We will prove that {g^l > {gsl on the border of two circles with centre in 
the origin and radii 1 and 4/5. 

In the following we first give a function G that is an upper bound for g, 
then we study separately the maximum of {gsl and the minimum of on 
the two circles. After that we prove the theorem. 

An upper bound for g. Let 

M := l-maxak,i = a2,i/2 

Z k,i 

(It follows from the fact that maxj ak^i = «fc,[fe/2] and that cek,[k/2] is decreasing 
in k.) 

We introduce another auxihary sequence {Dk}k- 
Do = do 

Di = di (20) 

Lemma 20 Let G{z) := '^'kLo^kZ^ ■ The radius of convergence of g is 
larger than the radius of convergence of G which is larger than [3"^ . 

Proof. Thanks to (fT4|) . by induction < d^ < for all k. 
By adding the third one of (l20ll we find 

k 

G{z)-do-diz = ^Dk+iz''+^ = Mz^^D,z'Dk-iZ^~' = MzG{zf-Mdlz 

k>l k>l i=0 
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hence for 2; e C inside the convergence radius of the series: 

1 ± ^1 - 4:Mz{do + {di - Mdl)z) 



G{z) 



2Mz 



The right sign is '— ' for aU z since G(0) = d^ and G is continuous. After 
some algebraic manipulations we write 

. ^ l-^/l-AMdoz + AAM'^dlz-' 

2Mz ^ ' 

where A = 1 - q;o,o/2M. 

The radius of convergence of G is the distance from the origin of the 
nearest zero of the radicand. The roots of the latter are both real and positive: 

The smallest one is the convergence radius of G and one can verify that 

(q;2,i - Q;o,o)"o 

This completes the proof. ■ 

On the convergence of g, g and the maximum of cib- 

Lemma 21 The radius of convergence of g and hence g is greater than 1. 
gB on the circumferences of radii 1 and 4/5 is bounded by: 

|^i?(e^'')| < 0.062, V^G[0,27r] 
|^s(4/5e*^) I < 0.031, e [0, 27r] 

Proof. Let z e C, 

00 00 
\~gB{z)\<2\z\J2K\\z\'^2j2K\ 

k=2 k=2 



Z 



\k+l 
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By the estimates on j, we have |1 — a^^il < CjS"^^; we also need some 
lower bound and, for A; > 1, j > ai q- Then for all A; > 1: 



l"fc+il 



\ Y^{ak,^ - l)didk-i < ]-C(3'^^ didk 



j=0 



i=0 



1 ^ 

< 7 ^ak,ididk- 



"1,0 2 



2=0 



"1,0 



"1,0 



Thanks to fl22l) . this proves that the radius of convergence of g is greater 
than 1. 

Putting together the last two upper bounds we get 



riff 

\~9Bi^)\<2y2— D,\z 



2k-2 



-1 Cf3- 



k=2 



"1,0 



"1,0 



k=2 



Cp-^\z\ 

"1,0 



{Gi/3^\z\) - Do - D,f3^\z\) 



The value of G{P'^\z\) can be computed thanks to (I2T1) . and the bounds one 
gets are those in the lemma. ■ 



On the minimum of gj^. 

Lemma 22 The maximum and the minimum of \gA{z)\ on the circumfer- 
ences with centre in the origin and radii 1 and 4/5 are on the following 
points: 



^^(-1) ^3.170 
^a(-4/5) ^ 2.650 



^^(1) ^ -0.092 
^a(4/5) ^ 0.040. 



Proof. We study \gA{pe'^)\- Recall that gA{z) = 1 - 2d'QZ - 2d[z^ 
Yl'k=o ^kz'^ is a polynomial with real coefficients. 



lUpe"')]' = ~gAipe'%Aipe^') = Upe^'TgA^pe-') = Yl ^^^^P^ 

j,k=0 



+k^ie{j-k) 



Y c]p^^ + 2(coCip + CiC2p^) cos(^) + 2coC2P^ cos(20) 
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By setting ^|(y(^(pe*^)p = we find: 

= (cqCi + ciC2p^) sin(0) + 2coC2psm{29) 
= (cqCi + CiC2P^ + AC0C2P cos{9)) sm{9) 

The first factor is never because |coCi + CiC2P^\ — |4coC2p| > 

|coCi + ciC2p^| - |4coC2p| = 2do|l - 2d[p'^\ - 8\d[\p 

= 2cio(l + (1 - «o,o)c^y ) - 4(1 - ao,o)dlp 
= 2do[l + (1 - ao,o)c?oP^ - 2(1 - ao,o)c?op] 
> 2cio[l + (1 - «o,o)^c?y - 2(1 - ao,o)dop] 

= 2rfo(l - (1 - ao,o)doP? > 

So the maximum and the minimum of |^a(-2^)| on the circumference with 
centre in the origin and radius p are on the real points ±p. 

Setting p = 1 and p = 4/5 and computing values we get the thesis. ■ 

Proof of Theorem 1191 Direct computation shows that qa has only one 
root inside the unit circle, in some real point in (4/5, 1). By Rouche's theorem 
and the bounds of lemmas [21] and [22] above, g has only one root inside the 
circular crown with centre in the origin and radii 4/5 and 1. This root is a 
real number R G (4/5, 1), since the same estimates state that (by continuity) 
g must be zero for some real number in the interval: 

^(1) = ^^(1) + ^b(I) < < ^a(4/5) + ^b(4/5) = ^(4/5). 
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